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ABSTRACT 

Aims. We aimed to monitor the optical linear polarimetric signal of the magnetized, rapidly rotating M8.5 dwarf TVLM 513-46546. 
Methods. R- and /-band linear polarimetry images were collected with the ALFOSC instrument of the 2.56-m Nordic Optical 
Telescope on two consecutive nights covering about 0.5 and 4 rotation cycles in the R and 1 filters, respectively. We also obtained 
simultaneous intensity curves by means of differential photometry. The typical precision of the data is ±0.46% (R), ±0.35% (/) in the 
linear polarization degree and ±9 mmag (//), ±1.6 mmag (/) in the differential intensity curves. 

Results. Strong and variable linear polarization is detected in the R and I filters, with values of maximum polarization (p* = 
1.30 ±0.35%) similar for both bands. The intensity and the polarimetric curves present a sinusoid-like pattern with a periodicity 
of ~1.98 h, which we ascribe to structures in TVLM513-46's surface synchronized with rotation. We found that the peaks of the in¬ 
tensity and polarimetric curves occur with a phase difference of 0.18 ± 0.01, and that the maximum of the linear polarization happens 
nearly half a period (0.59 ± 0.03) after the radio pulse. We discussed different scenarios to account for the observed properties of the 
light curves. 
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1. Introduction 


The M8.5 dwarf TVLM 513-46546 (TVLM 513-46: Tinney 
et al. 1993!) is one of the most studied ultracool dwarfs in 
the literature. With a trig onometric distance of 10.76 ± 0.03 
pc IlLorbrich et alJ 20131), an effective temperature (T e g) of 
2175 ± 150 K, a l umino sity of log (L/L©) = -3.59 + 0.02 dex 
(iGolimow ski et all 20041), and a lack of Li at 670.8 nm ( Martin 


et al. l 19941 : Reid et al.! 2002) the models of lBaraffe et all (20031) 

constrain the mass of TVLM 513—46 to the range 0.06-0.08 Mg 
and its radius to 0.1 R© for ages older than 0.5 Gyr. It shows Her 
in emission that changes moderately with time, which suggests 
some chromospheric activity (IReid et al.l 20021) Lrom its high 
projectedjxffation^ektcit^vsin^^JjO^^JcmvA^Mohanty 
& Basri 20031) and its rota tion period (1.959574 ± 0.000002 h, 
IWolszczan & Routel 20141), the inclination of its rotation axis is 
found to be i = 74.5^° g 2 deg. The detection of radio emission 
reveal s a multipolar magnetic field, with inten sities as high as 
3 kG (lHallinan et al.l 120061 : IBerger et alJ 20081) This magnetic 
field and the detection of optical variability have led to a rich 
discussion on which physical scenario is the most likely to ex 


plain its variability: magnetic spots or dust clou ds (Lane et al 


2007; Littlefai r et alJ 120081 : IBerger et al.l 120081 : lHarding et al 


2013D 


Because the identification of this scenario is degenerate, 
additional observational parameters which add new informa¬ 
tion from this objects’s atmosphere are necessary. In particu¬ 
lar, the amount of linear polarization suffered by the light emit¬ 
ted/scattered toward a distant observer, and its variability, could 


be a very powerful discrimination tool. In this work, we present 
our R- and /-band monitoring of the intensity and linear polari¬ 
metric properties for TVLM 513—46 to try to shed more light on 
the physical processes responsible for the optical modulation. 


2. Observations 

We collected R- (/l c =631 nm) and /-band (A r = 810 nm) lin¬ 
ear polarimetry images of TVLM 513-46 using the Andalucla 
Laint Object Spectrograph and Camera (ALLOSC) mounted on 
the 2.56-m Nordic Optical Telescope (NOT) on 2013 May 18 
and 19. The target was monitored over ~3.5 and ~4.2 h (/-band) 
during the first and second nights, and over ~1 h (//-band) during 
the second night. We thus covered ~ 4 rotation cycles in the I fil¬ 
ter and half a rotation cycle in the //-band with a typical cadence 
of 6.2 and 9.6 min. Using these data, we also retrieved the R- and 
/-band intensity curves of TVLM 513-46 by means of differen¬ 
tial photometry and by combining the polarimetric ordinary and 
extraordinary rays. The star 2MASS J15011008+2250069 (2.05 
mag brighter in I and located near our target) acted as the ref¬ 
erence/comparison source. The detailed description of the ob¬ 
serving strategy, instrumental configuration, and data reduction 
is provided in Sect. [A] (online material). The normalized Stokes 
parameters q and n, the degree of linear polarization (/*), and the 
polarization vibration angle (0) were derived using the flux-ratio 
method and equations 1-4 from [Zapatero Osorio et al] (2005D 
The final data have typical uncertainties as follows: ±0.46% (//) 
and +0.35% (/) in P, +9° in 0, and +3.4 mmag (/, first night), 
+ 1.6 mmag (/, second night) and +9 mmag (R) in the differential 
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Fig. 1. Stokes q - u plane for the /-band measurements of 
TVLM513-46 (red circles) taken during the second night of 
observations. Instrumental polarization is shown with a white 
circle. The typical uncertainties in q and u are represented in the 
top-left. 



Fig. 2. Differential intensity (top) and debiased linear polarime- 
try (bottom) curves of TVLM513-46 taken on 2013 May 18 
(left) and 19 (right). The vertical dashed lines separate the I- 
(blue circles) and the /Lband data (gray circles). TVLM 513-46 
is shown with solid dots. The reference star is depicted with 
white triangles and no error bars (for clarity) in the bottom panel. 
The horizontal dashed line stands for the upper limit on instru¬ 
mental polarization. Zero time corresponds to Julian dates of 
2456431.5734 (2013 May 18) and 2456432.4922 (2013 May 
19). 


intensity curves. Figure Q] illustrates the /-band q - it plane ob¬ 
tained using the second observing night measurements. We shall 
focus on these data because they have better quality by a fac¬ 
tor of 2 than the photometry of the first night (see also Fig. [2]). 
This allows us to study the details of the polarimetric light curve 
of TVLM 513-46. The data of the first observing night are pre¬ 
sented in Fig.[2]and Section A; they support our findings within 
the quoted error bars. 


3. Light curves 

Figure [2] (top panel) illustrates the R- and /-band differential 
intensity curves of TVLM 513—46 as a function of time. The 


photometry exhibits sinusoidal va riabi lity, a Lomb-Scargle peri- 
odogram of our curved JScarglel 198211 shows a peak at 1.9 + 0.4 
h (Figure IA.3I ) , which fully agrees with the period reported by 
other groups ( Wolszczan & Routell2014 . and references therein) 
and confirms the stable periodic variability status of TVLM 
513—46. This period likely reflects the rotational modulation of 
TVLM 513-46. The error bar of +0.4 h was determined as the 
full width at half maximum of the 1,9-h peak of the periodogram. 
The /-band amplitudes of the differential intensity curves are 
measured at ~4.9 mmag (first night) and ~3.9 mmag (second 
night), which are interm ediate between tho se gi ven for filters of 
related wavelengths by lLane et akl (l2007t) and iLittlefair et alJ 
(2008!) Our measurements support the changing amplitudes of 
the observed sinusoidal variations with time. Because the peri¬ 
odicity is constant, this suggests that the intensity and/or size of 
the feature responsible for the variations is also changing (see 
also IHarding et alJl2013l) . 

The bottom panel of Fig. [2] depicts the evolution of the R- 
and I -band linear polarimetry of TVLM 513-46. Because of the 
better signal-to-noise ratio of the second night photometry, the 
modulation of the polarimetric light curve becomes apparent on 
the right panel. Values between zero and strong linear polariza¬ 
tion of p* » 1.30% + 0.35% are seen at the two wavelengths 
of our study. As a consistency check, we also derived the linear 
polarimetry photometry of the reference star in the same man¬ 
ner as for the target. These measurements are included in Fig. [2] 
(bottom panel); they uniformly lie below the upper limit on the 
instrumental polarization, as expected for an unpolarized source. 
Because the reference star is located at a projected separation of 
~Y1" from our dwarf, we concluded that the observed high lin¬ 
ear polarization is intrinsic to TVLM 513—46 and is not due to 
an interloper cloud. Furthermore, the dwarf lies at a distance of 
10.76 ± 0.03 pc ( Forbrich et al J2013h in a direction with negligi¬ 
ble extinction by the interstellar medium. The solar vicinity (<25 
pc) shows very little linear polarization (less than 0.0004% per 
parsec, Piirola 1977). Only at dist ances >70 pc, inters tellar po¬ 
larization becomes non-negligible (iTamburini et al.1 2002!) 

Interestingly, the maximum polarization signal of TVLM 
513—46 is similar for both the R and / filters, although we cau¬ 
tion that the //-band observations are not simultaneous with the 
/-band and that the former data cover only a fraction of the ro¬ 
tation period. The Lomb-Scargle periodogram of the /-band lin¬ 
ear polarimetry light curve shows a peak at ~2-h (Figure IA.3b 
with a confidence of 83% (see the Appendix). This agrees with 
the periodicity previously derived from the intensity curve and 
with t he v ery precise rotation period s given by IHarding et alJ 
(1201 3l) and lWolszczan & Routel (120141) . A similar Lomb-Scargle 
analysis of the polarimetric photometry of the reference star 
yielded no significant peaks at around 2 h. Our data indicate 
that TVLM 513-46 is linearly polarized and that the intensity 
of the polarization changes with a periodicity compatible with 
the dwarf’s rotation. 

The intensity and linear polarimetry light curves of the M8.5 
dwarf are folded in phase using a period of 1.9798 h in Fig. [3] 
(and in Fig. IA.4l for data of both nights). Despite the short-time 
coverage of the //-band observations and their relatively large 
associated uncertainty, our data indicate that the intensity vari¬ 
ations in the R and / wavelengths are correlated (Fig. IA.4L as 
opposed to the anti-phased g and i light curves discussed by 
ILittlefair et al.1 (2008!) The previous 7-band (1.2 pm) linear po- 


1 We used our own codes and the codes provided by the NASA 
exoplanet archive on: http://exoplanetarchive.ipac.caltech.edu/cgi- 
bin/Periodogram/nph- simpleupload 
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Phase 

Fig. 3. /-band differential intensity (top), linear polarimetry 
(middle), and polarization vibration angle (bottom) curves 
folded in phase using a periodicity of 1.9798 h. In the top panel, 
each data point stands for the average of four individual mea¬ 
surements and its associated dispersion. In the bottom panel, 0 
is plotted for P/crp > 2.7. Two phases are presented for clarity. 
The green triangle stand s for a 7-band measurement taken from 
Miles-Pae z et al.1 (201 31) 


larimetry measurement by our group dMiles-Paez et al.l 20131) is 
also included. It was taken 3.7 months prior to the ALFOSC ob¬ 
servations. To improve the S/N ratio of the I -band intensity data 
and to equalize the time sampling of both light curves, we plot 
the average of four individual data points in the top panel. Within 
the observed scatter, the intensity and polarimetric light curves 
exhibit sinusoid-like variations, which we ascribe to large-scale 
structures in the dwarf’s surface that move with rotation. There 
are two striking properties of the linear polarimetry light curve: 
one is that the maximum of the I -band linear polarization oc¬ 
curs after the peak of the intensity curve. By applying simple 
sinusoidal fits to the data, we derived that both peaks are phase 
shifted by 0.18 + 0.01 or about 21 min. Wols zczan & Routel 
(20141) reported that the radio emission bursts of TVLM 513-46 
are detected 0.41 ± 0.02 of the period earlier than the peak of the 
intensity curve, which indicates the dipolar nature of the mag¬ 
netic field and the inclination of thejmagnetic field with respect 
to the rotation axis dBerger et al.l 2009) Therefore, there is a de¬ 
lay of nearly half a rotation (0.59 + 0.03) between one burst in 
the radio emission of the M8.5 dwarf and its /-band polarization 
maximum. 

The highest /-band linear polarization degrees take place 
during phases ~0.4 through ~0.8, while the intensity light curve 
keeps decreasing toward the minimum brightness (top and mid¬ 
dle panels of Fig. [3}. The polarization vibration angle is shown 
folded in phase in the bottom panel of Fig. [3 With few ex¬ 
ceptions, all I -band measurements are consistent with the mean 
value of 0 = 31°±10°, thus we do not resolve any particular 
pattern in the distribution of © given the precision of the data. 
Furthermore, the I- and 7-band polarization angles are compati¬ 
ble at the 1-cr with the quoted uncertainties, which suggests the 


little dependence of 0 on wavelength and/or the stability of the 
polarizing feature for at least 3.7 months. 


4. Interpretation 


The shape and characteristics of the polarimetric and intensity 
light curves of TVLM 513 -46, combined with other well-known 
properties of the M8.5 dwarf, may provide evidence of the phys¬ 
ical origins of its linear polarization. Because of the coincidental 
periodicities between the two data sets, it is likely that the po¬ 
larimetric variability is related to structures that spin with the 
dwarf’s rotation. For example, a (cloudy or magnetic) spot on 
the dwarf’s surface can change the intensity curve in a periodic 
manner, which will also be reflected in the polarization. Next, we 
discuss the following linear polarization mechanisms: Magnetic 
Zeeman splitting, scattering processes by electrons and dusty 
particles, and dichroic extinction. 

TVLM 513-46 harbors a large scale, steady magnetic field, 
which is responsible for the observed radio emission (Hallinan 
et al. 120061 : i Bergeret alJ20081) iMenardet al. (200 2) and Menard 
& Delfosse (20041) studied 21 M1-L5 dwarfs located at distances 
of less than 32 pc, finding values of P < 0.2% (R and / fil¬ 
ters). We remark that some of the stars analyzed by these au- 
thors have magnetic fields with strengths of a few to several kG 
dDonati et al.l2006i : lReiners & Basril 2007) similar in intensity to 
TVLM513-46. Nevertheless, magnetic Zeeman broadening and 
splitting of individual lines due to species like FeH and CrH in 
very cool M and L dwarfs have been observed and modeled theo¬ 
retically fe.g jReiners & Basri 2006; Kuzm vchov & Berdv ugina 
120131) . This yields significant linear polarization as a function of 
the magnetic field strength at very specific wavelengths. When 
convolved with high rotation velocities and broad-band filters, 
the polarimetric signal is notoriously reduced but the net integra¬ 
tion over broad bands may be different from zero, as pointed out 
by [Kuzmvchov & Berdvugina! (20131) The /-band filter of our 
observations covers features due to Nai, TiO, VO, CrH, FeH, 
HjO, ... in the wavelength interval 0.720-0.875pm, some of 
which are magnetically sensitive; therefore, we cannot discard 
a small contribution of this effect to the observed polarization. 
However, it may contribute with less than P 0.2% as sug¬ 
gested by the observations of ^emntf^nfljjyOO^ and Menard 
& Delfosse (2004!) As we detect very high values of linear po¬ 
larization spanning over 40% of the rotation period, the Zeeman 
broadening/splitting due to a strong magnetic field is not suffi¬ 
cient to account for all of the observed optical linear polarization 
of TVLM 513-46. 


In addition, the distribution of our polarimetric measure¬ 
ments in the q-u plane of Fig.Qjhints at polarizing mechanisms 
rel ated to scattering p r ocesse s. Such a distribution is reproduced 
bv IClarke & McGala (19861) who modeled the polarimetric be¬ 
havior of a globule of particles in orbit around a central unpo¬ 
larized source using the Thompson scattering scenario (see their 
Figs. 3b and 5a). For an equator-on source and a cloud of parti¬ 
cles located close to this region, the Stokes parameters describe 
a near circular pattern like that of Fig. [T|( also Fig. lA.2l >. The q-u 
distribution shape changes significantly as a function of the lo¬ 
cation of the cloud across the disk of the source and the angle 
between the magnetic and rotational poles. The strong magnetic 
field of TVLM513-46 reveals the presence of a surrounding 
envelope of low pressure electrons. Actually, the 100% circu¬ 
larly polarized radio pulses of the M8.5 dwarf are believed to 
be most likely caused by the electron-cyclotron maser instabil¬ 
ity (lHallinan et al,ll2008l ), while the quiescent radio emission is 
probably due to gyrosynchrotron radiation. Our data are consis- 
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tent with the R- and /-band linear polarization being originated 
by the Thompson scattering of an assembly of free electrons 
in a similar manne r to what is observed in h ot Be stars (e.g., 
iBrown et al.l Il978k IClarke & McGaldll986l) . Under this sce¬ 
nario, the electronic cloud must be located close to one radius 
of TVLM 513-46 t o produce polari metri c variability synchro- 
nized with rotation. lYu et al.l ( 2011 1 and iKuznetsov & Vlasov! 
(12012h determined that the electron density is of the order of 
10 5 cirT 3 in TVLM513-46’s magnetosphere, which is confined 
within two stellar radii according to lJaeger et al.l (2011!) 

Dust, which is present in primordial and debris disks, can 
also act as an efficient polarizer. Infrared flux excesses provide 
evidence of the presence of any of these disks surrounding cen¬ 
tral objects. TVLM 513-46 does not appear to have clear flux 
excesses at 24 pm dGautier et al.l 20071), which thus discards the 
existence of warm primordial and debris disks. In addition, the 
M8.5 dwarf has efficiently depleted lithium from its atmosphere 
dReid et alJ 2002T), which indicates a minim um age of a few hun¬ 
dred Myr (also see lHallinan et al.1l2008ll . Dust is also present 
in the phot ospher es of dwarfs cooler than T e g = 2100 K (e.g., 
iTsuii et ail 1996!l as a result of the natural chemistry of low- 
temperature and high-gravity atmospheres that includes the con¬ 
densation of refractory elements into liquid and solid particles. 
Under the scenario of large rotation velocities (that induce oblate 
shapes) and/or inho mogeneous distributions of dusty cl ouds, 
ISengupta & Marlevl (201 01) and lMarlev & Sengupfal (201 1|) the¬ 
oretically showed that the net linear polarization of dusty dwarfs 
by single scattering processes can be as hig h as a few per cent at 
optical and near-infrared wavelengths. iMiles-Paez et al.l (120131 
and references therein) provided observational proof of linear 
polarization detection s in ultracool dwarfs. TVLM 513-46 has 
T e ff = 2025-2480 K dGolimowski et all l2004t Isteohens et alJ 
2009) and a rapid rotation (1.96 h and v sin i = 60 kms '); con¬ 
densates (dust particles) are thus likely populating the upper 
(and coolest) atmospheric layers. If these particles are unevenly 
distributed across the atmosphere and/or aligned with the in¬ 
tense magnetic field, some linear polarimetry signal would be 
expected. 

Ide Kok et alJ (201 ID modeled the near-infrared polarization 
of a rotating dwarf with a 20°x20° dusty spot on its equator 
(Rayleigh scattering). They found a varying polarization degree 
and vibration angle with amplitudes of ~0.6% and ~40° (A - 
1.11 pm) as the cloudy spot moves across the object’s visible 
disk with resgectJoJhe_observer 1- A^_illustrated > hi - Fi i g il _5 > b^ i de 
Kok et al. (201 ID, the polarization degree has two maxima sep¬ 
arated by a phase difference of 0.17 and a deep minimum in- 
between just when the cloud lies along the line of sight. There is 
a hint of this feature at phase 0.5 in the middle panel of Fig. [3 
Regarding 0, the model predicts a sudden change in the angle 
within a phase difference of 0.055, which translates into 6.6 min 
in the case of the M8.5 dwarf. Our data do not have enough 
temp oral sampling to test these nredictions in p* and Q. dc Kok 
et al. (201 ID also calculated the intensity curve caused by the 
equatorial dusty spot. It shows a mild peak at the same phase 
as the deep minimum between the two polarimetric maxima. 
Our TVLM 513-46 data do not recreate this signature since the 
maximum brightness of the /-band differential light curve and 
the deep minimum between the two maxima of the polariza¬ 
tion curve are separated by ~0.15 of the rotation period. This 
phase difference might be explained by changing the geometry 
in which the grains are illuminated. 

Finally, dichroic extinction by magnetically aligned non- 
spherical grains ma y lead to significant pola rization ( Davis & 
Greenstein 195 ID; this mechanism is believed to cause the inter¬ 


stellar polarization observed toward many stars in our Galaxy 
and has_been used to map out the Galactic magnetic field 
(iMathewson & FordT 197(1) Under this scenario, the linear polar¬ 
ization of TVLM 513-46 would lie parallel to the magnetic field 
direction, which is 31 + 10°. The periodic and sinusoid-like vari¬ 
ation of p* would be caused by a cloud of aligned non-spherical 
grains rotating with the dwarf. The impact of convective motions 
and rapid rotation on such dusty clouds is not clear. Models of 
dichroic polarization particularized for ultra-cool dwarfs are re¬ 
quired to explore this scenario further. 

From our discussion it is clear that, with the present data, it 
is not posible to solve unambiguously the origin of the observed 
optical linear polarization. To this end, follow-up observations 
taken at shorter and longer wavelengths and with a higher ca¬ 
dence would also be beneficial: On the one hand, the intensity 
of the dichroic polarization has a strong dependence on wave- 
length, whereas the polariza tion vibration angle keeps constant 
dTreanoill 1963b ICovnd 1974D , and on the other hand, linear po¬ 
larization due to other mechanisms shows a moderate variation 
with wavelength. Additionally, polarimetric observations com¬ 
bined with other techniques to map out the presence of dusty 
clouds, such as the tomography of the surfaces (ICrossfield et al.l 
2014D could be valuable tools to determine the geometry, distri¬ 
bution, and other properties of the clouds in the atmospheres of 
ultracool dwarfs, especially for those that do not harbor strong 
magnetic fields such as TVLM513-46. 

Acknowledgements. This work was based on observations made with the Nordic 
Optical Telescope operated on the island of La Palma by the Nordic Optical 
Telescope Scientific Association (NOTSA) in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astroffsica de Canarias. This work 
is partly financed by the Spanish Ministry of Economics and Competitiveness 
through projects AYA2012-39612-C03-02 and AYA2011-30147-C03-03. 

References 

Baraffe, I., Chabrier, G.. Barman, T. S., et al. 2003, A&A, 402, 701 
Berger, E., Gizis, J. E., Giampapa, M. S., et al. 2008, ApJ, 673, 1080 
Berger, E., Rutledge, R. E., Phan-Bao, N., et al. 2009, ApJ, 695, 310 
Brown, J. C., McLean, I. S., & Emslie, A. G. 1978, A&A, 68, 415 
Clarke, D. & McGale, P. A. 1986, A&A, 169, 251 
Coyne, G. V. 1974, AJ, 79, 565 

Crossfield, I. J. M., Biller, B., Schlieder, J. E., et al. 2014, Nature, 505, 654 
Davis, Jr., L. & Greenstein, J. L. 1951, ApJ, 114, 206 
de Kok, R. J., Starn, D. M„ & Karalidi, T. 2011, ApJ, 741, 59 
Deeg, H. J. & Doyle, L. R. 2001, in TWP, ed. W. J. Borucki & L. E. Lasher, 85 
Donati, J.-F., Forveille, T., Collier Cameron, A., et al. 2006, Science, 311, 633 
Forbrich, J., Berger, E., & Reid, M. J. 2013, ApJ, 777, 70 
Fossati, L., Bagnulo, S.. Mason, E., et al. 2007, in ASPCS, Vol. 364, , 503 
Gautier, ID, T. N., Rieke, G. H., Stansberry, J., et al. 2007, ApJ, 667, 527 
Golimowski, D. A., Leggett, S. K., Marley, M. S., et al. 2004, AJ, 127, 3516 
Hallinan, G., Antonova, A., Doyle, J. G., et al. 2006, ApJ, 653, 690 
Hallinan, G., Antonova, A., Doyle, J. G., et al. 2008, ApJ, 684, 644 
Harding, L. K„ Hallinan, G., Boyle, R. P„ et al. 2013, ApJ, 779, 101 
Home, J. H. & Baliunas, S. L. 1986, ApJ, 302, 757 
Jaeger, T. R., Osten, R. A., Lazio, T. J., et al. 2011. AJ. 142, 189 
Kuzmychov, O. & Berdyugina, S. V. 2013, A&A, 558, ADO 
Kuznetsov, A. A. & Vlasov, V. G. 2012, A&A, 539, A141 
Lane, C., Hallinan, G., Zavala, R. T., et al. 2007, ApJ, 668, LI63 
Littlefair, S. P., Dhillon, V. S., Marsh, T. R„ et al. 2008, MNRAS, 391, L88 
Lomb. N. R. 1976, Ap&SS, 39, 447 
Marley, M. S. & Sengupta, S. 2011, MNRAS, 417, 2874 
Martin, E. L., Rebolo, R., & Magazzu, A. 1994, ApJ, 436, 262 
Mathewson, D. S. & Ford, V. L. 1970, MmRAS, 74, 139 
Menard, F. & Delfosse, X. 2004, in SF2A-2004, 305 
Menard, F.. Delfosse, X., & Monin, J.-L. 2002, A&A, 396, L35 
Miles-Paez, P. A., Zapatero Osorio, M. R.. Palle, E., et al. 2013, A&A, 556, 
A125 

Mohanty, S. & Basri, G. 2003, ApJ, 583, 451 

Mortier, A., Faria, J. P. Correia, C. M., et al. 2015, A&A, 573, A101 

Piirola, V. 1977, A&AS, 30, 213 


4 
















































P. A. Miles-Paez et al.: Rotational modulation of the linear polarimetric variability of the cool dwarf TVLM 513—46546 

Reid, I. N„ Kirkpatrick, J. D„ Liebert, J„ et al. 2002, AJ, 124, 519 

Reiners, A. & Basri, G. 2006, ApJ, 644, 497 

Reiners, A. & Basri, G. 2007, ApJ, 656, 1121 

Scargle, J. D. 1982, ApJ, 263, 835 

Sengupta, S. & Marley, M. S. 2010, ApJ, 722, L142 

Stephens, D. C., Leggett, S. K., Cushing, M. C., et al. 2009, ApJ, 702, 154 

Tamburini, F„ Ortolani, S., & Bianchini, A. 2002, A&A, 394, 675 

Tinney, C. G„ Mould, J. R„ & Reid, I. N. 1993, AJ, 105, 1045 

Treanor, P. J. 1963, AJ, 68, 185 

Tsuji, T„ Ohnaka, K„ & Aoki, W. 1996, A&A, 305, LI 

Wardle, J. F. C. & Kronberg, P. P. 1974, ApJ, 194, 249 

Whittet, D. C. B„ Martin, P. G„ Hough, J. H„ et al. 1992, ApJ, 386, 562 

Wolszczan, A. & Route, M. 2014, ApJ, 788, 23 

Yu, S„ Hallinan, G„ Doyle, J. G„ et al. 2011, A&A, 525, A39 

Zapatero Osorio, M. R., Caballero, J. A., & Bejar, V. J. S. 2005, ApJ, 621, 445 


5 



P. A. Miles-Paez et al.: Rotational modulation of the linear polarimetric variability of the cool dwarf TVLM 513^16546, Online Material p 1 


Appendix A: Observations and data reduction 


We collected linear polarimetry images of TVLM 513-46 
using the Andalucla Faint Object Spectrograph and Camera 
(ALFOSC) mounted on the 2.56-m Nordic Optical Telescope 
(NOT) on 2013 May 18 and 19. ALFOSC has a 2048x2048 
E2V detector with a pixel size of O'/19. The target was moni¬ 
tored in the /-band (A c = 810 nm) over ~ 3.5 and ~ 4.2 h during 
the first and second nights, and in the //-band (A c = 631 nm) 
over ~ 1 h during the second night. We thus covered ~ 4 ro¬ 
tation cycles in the /-band and half a cycle in the //-band. Sky 
conditions were clear and seeing varied in the interval T/l-l"5. 
Observations were carried out at air masses ranging from 1.05 
through 2.34. 

The linear polarimetry observing mode of ALFOSC consists 
of a half-wave plate and a calcite block, which provides simulta¬ 
neous images of the ordinary and the extraordinary beams sep¬ 
arated by 15". The total unvignetted field of view is 140" in di¬ 
ameter, which we rotated by 103° East of North to align TVLM 
513-46 and a nearby, bright reference star along the v-axis of 
the detector, as it is indicated in Fig. IA. II Individual exposure 
times were 50 s (/, first night) and 100 s (R and /, second night) 
per position of the half wave plate (0°, 22/5, 45°, and 67/5). On 
May 18, polarimetric images were taken at two nod positions on 
the detector separated vertically by 10" for a proper sky subtrac¬ 
tion. This nodding pattern was not applied on the following night 
since TVLM 513-46 is detected well above the sky contribution 
with signal-to-noise (S/N) ratio of ~220 and ~300 (/-band) in 
50-s and 100-s individual integrations, respectively (S/N is mea¬ 
sured at the peak of the energy distribution of the source im¬ 
ages). One polarimetric cycle was completed every 6.2 min (/, 
first night) and 9.6 min (R and /, second night), including over¬ 
heads. This allowed us to sample one rotation of TVLM 513-46 
using a minimum of 12 and a maximum of 19 linear polarimetry 
measurements. Raw images were bias subtracted and flat-field 
corrected using data acquired during dawn and dusk and routines 
within IRAffl. We also observed non-polarized stars (GJ 838.4 
and WD 1615-15 4.lFossati et al.ll2007t) and one polarized source 
(Hilter 997. IWhittet et al. 11992 ) by employing the same instru¬ 
mental configuration and on the same observing dates as the tar¬ 
get. This allowed us to check for the efficiency of the instrument 
and to set an upper limit on the instrumental polarization. 

We derived the normalized Stokes parameters q and u, the 
degree of linear polarization (P), and the polarization vibration 
angl e (0) using the flux-rat io method and equations 1-4 from 
IZapatero Osorio et al.l (20051) The Stokes parameters and the an¬ 
gle 0 were properly corrected for the rotation of the field of view. 
Fluxes of each ordinary and extraordinary images were extracted 
by defining circular apertures ranging between 0.5 and 6 times 
the size of the full-width-at-half-maximum (FWHM) and using 
the PHOT package within IRAF. Sky annuli had inner radii rang¬ 
ing from 3.5 x to 6 xFWHM and widths of 1.5 xF WHM . The 
detailed procedure is described in Miles-Pa ez et al.l (20131) Our 
final q, u, P, and © values result from averaging the aperture pho¬ 
tometry in the interval 2-4 x FWHM; the errors associated with 
q and u are defined as the standard deviations of the measure¬ 
ments in the selected range of apertures, the errors in P are com¬ 
puted as the quadratic sum of the q and u quoted uncertainties 
plus the uncertainty introduced by a possible instrumental linear 
polarization (see below), and errors in © are determined follow- 


1 IRAF (Image Reduction and Analysis Facility) is distributed by the 
National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co¬ 
operative agreement with the National Science Foundation. 



Fig.A.l. ALFOSC linear polarimetric I -band image of 
TVLM 513-46. The circle indicates the unvignetted area. 
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Fig.A.2. Stokes q-u plane for the R- and /-band measurements 
of TVLM 513-46) taken during first and second night of obser¬ 
vations. Instrumental polarization is shown with a white circle. 
The typical uncertainties in q and u are indicated. 


ing the equations given in Miles^Paez^etjd. (2013) and Wardle 
& Kronberg (19741) We caution that the expression for the 0 
uncertainty is valid for P/cr P > 3, where cr P is the error in lin¬ 
ear polarimetry. Our linear polarimetric data of TVLM513-46 
have a typical uncertainty of +0.46% (R) and ±0.35% (/) in P, 
and ±9° in 0. 
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Fig. A.3. Top: Lomb-Scargle periodograms of the intensity (left) 
and polarimetric (right) light curves of TVLM 513-46 (red) and 
the reference star (black). Bottom: Relative probabilities of the 
explored periods using a Bayesian formalism for the generalized 
Lomb-Scargle periodogram (see text). 


Using the observations of the zero-polarized standard stars, 
we checked that the instrumental polarization lies below 
0.31% and 0.27% in the R- and /-bands. The polarized stan¬ 
dard star allowed us to determine the zero point correction 
for the polarization vibration angle to be 0 o ,r = (F8+2'/8 and 
0o./ = L3+3/0, which agrees with the values tabulated in the 
ALFOSC manuaQ. From now on, we shall use the debiased lin¬ 
ear polarization degree, p* , defined as 

P* = ^ P 2 -cr 2 p , (A.l) 

which accounts for an overestimation of the polarizationjdgnal 
at low S/N or small values of P/crp (tWardle & Kronbergi 19741) 
We set p* = P if cr P > P. In Figs.[2][3] and lA.4l these polarimet¬ 
ric data are displayed with an associated arrow indicating that 
the true polarization index likely lies between 0% and the sym¬ 
bol position. We note that the dispersion of consecutive © mea¬ 
surements is of the order of the uncertainty in the polarization 
vibration angle when P/crp > 2.7, which is close to the regime 
where th e calculation of 0 has statistical significance ( Wardle & 
Kronberg 19741) This provides support to our method for deriv¬ 
ing the error bars associated with q, u, and P or p*. The individ¬ 
ual Stokes parameters values for the R and / filters are shown in 
the q - u plane of Fig. 1701 

In addition to the linear polarimetric evolution, we retrieved 
the R- and /-band intensity curves of TVLM 513—46 by means of 
differential photometry. The star 2MASS J15011008+2250069 
(2.05 mag brighter in / and located near our target) acted as 
the reference/comparison source. This is the only bright star 
in the ALFOSC field of view that was usable for the differen¬ 
tial photometry technique (Fig. A.l). Because of its lower lu¬ 
minosity, the error bars in the differential intensity curves are 
dominated by the photon noise of TVLM 513—46. The typically 
1.0 x FWHM-aperture fluxes of the ordinary and extraordinary 
images of TVLM 513—46 and its_ reference star were extracted 
with the VAPHOT package dDeeg & Dovlel 2001l). which works 
under the IRAF environment and is optimized for differential 

http://www.not.iac.es/instruments/alfosc/polarimetry/index.html 



Fig. A.4. R- and /-band differential intensity (top), linear po- 
larimetry (middle), and polarization vibration angle (bottom) 
curves folded in phase using a periodicity of 1.9798 h. Symbols 
as in Fig. IA.2I In the top panel, each data point stands for the 
average of four individual measurements and its associated dis¬ 
persion. In the bottom panel, © is plotted for P/crp > 2.7. Two 
phases are presented for clarity. 


light curves. To build the intensity curves, ordinary and extraor¬ 
dinary fluxes were combined per individual frame, which thus 
provides a time sampling four times higher than the cadence of 
the linear polarimetry data. We estimated the uncertainties in the 
differential light curves to be ±3.4 mmag (/-band) for the first- 
night photometry, and +1.6 mmag (/) and +9 mmag (R) for the 
second night data. We checked that the measured linear polar¬ 
ization degree and the differential intensity light curves do not 
correlate with airmass. 

The Lomb-Scargle periodograms of the /-band intensity and 
polarimetric light curves of T VLM 513 -46 computed following 
iLombl ( 19761) and lScarglel ( 19821) are illustrated in the top panels 
of Figure IA.3I The periodogram of the polarimetric light curve 
of the reference star is also included in the Figure. We explored 
frequencies between io = 2n/T and to = 2 N a /T with a spacing 
of about 1/T dHorne & Baliunas] 19861), where N a is the number 
of data points of the light curves (296 for the intensity curve, and 
74 for the polarimetric curve), and T is the total time coverage 
of the data (26.4 h). We thus surveyed periods in the range 0.3- 
26.4 h using 433 and 88 independent frequencies for the inten¬ 
sity and polarimetric data, respectively. The confidence levels of 
the 2-h peaks were estimated using the false-alarm-probability 
(1-FAP) at 99% (intensity light curve) and 83% (polarimetry). 
The bottom panels of Figure lA3l denict the relative probabilities 
according to the Bayesian formalism described in lMortier et al.l 
(20151) 
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